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We have carried out simulations of the hadronization 
of a hot, ideal but effectively massive quark-gluon gas 
into color neutral clusters in the framework of the semi- 
classical SU(3) chromodielectric model. We have studied 
the possible quark-gluon compositions of clusters as well 
as the dual mass distribution and spectra, aiming to ob¬ 
tain an insight into relations between hadronic spectral 
properties and the confinement mechanism in this model, 
pacs: 12.38.Mh, 12.38.Lg, 24.85-|-p 


I. INTRODUCTION 

The quest for the quark-gluon plasma is one of the ma¬ 
jor scientific efforts of the last decade of strong interaction 
physics. In recent years 200 AGeV Pb beams became 
available for the experimental study at CERN SPS, and 
future colliders will reach even higher energies, namely 
200 -b 200 A GeV at RHIG Brookhaven and 900 -b 900 
AGeV at GERN LHG. These energies correspond to 
an estimated maximal energy density of e = 5 GeV/fm^ 
and e = 9 GeV/fm^ or temperature T = 230 MeV and 
T = 260 MeV, respectively. At GERN SPS the achieved 
maximal temperature is about T = 160 MeV Q. 

Heavy ions like Au and Pb are the best means for ap¬ 
proaching the thermodynamical limit in particle accel¬ 
erator experiments as closely as possible. They provide 
high energy- and baryon densities. These advantages on 
the experimental side lead to trouble for theoretical stud¬ 
ies: although the thermodynamical approach to quark 
matter is highly developed a dynamical insight or 

finite baryon density studies are at their very beginning 
The solution of both theoretical problems, i.e. the 
real-time description and the study at finite density, is 
connected to the invention of new non-perturbative tech¬ 
niques, which are able to handle complex effective ac¬ 
tions. On the other hand, there are several phenomeno¬ 
logical models aiming to simulate these unknown dynam¬ 
ical features of quantum chromodynamics (QGD). 


^part of the Ph. D. thesis of C. T. Traxler 


In theoretical works the main goal is to understand 
the mechanism of color confinement. In this perspec¬ 
tive, understanding hadronization means learning about 
the confinement mechanism in an experimentally control¬ 
lable situation. We briefly review hadronization mod¬ 
els presently known to us grouped in three categories: 
prompt, equilibrium and evolutionary. The process takes 
zero (i.e. less than 1 fm/c), infinite (more than 20 fm/c) 
or a finite amount of time (2-6 fm/c) in these three 
respective categories. 

String- and parton-fragmentation models in the 
first category are the oldest prompt hadronization mod¬ 
els; they have been used already in the description of 
elementary particle collisions (pp, pp and e“'"e“). With 
the help of Glauber’s model g the yields of elementary 
pp-events are then summed up for a complicated heavy 
ion collision. In these calculations hadron ratios depend 
on phenomenoli^ical fragmentation ratios or fragmenta¬ 
tion functions M. Although these parameters are non- 
perturbative, any explanation or derivation of them from 
QGD is missing. 

Further statistical and combinatorical models connect 
quark numbers and (pre-)hadron numbers by assuming 
an underlying coalescence. The quark cluster model used 
e.g. in the parton cascade model or partially in 
string models like VENUS and (UjRQMD is sen¬ 
sitive to quark properties, like effective masses. The AL¬ 
GOR model assumes an extremely fast and far-from- 
equilibrium hadronization of massive quarks after gluons 
are fragmented. Although individual hadron formation 
is immediate, the formation time instant may vary over 
a range of several tens of fm/c in these approaches. In 
these models the hadron formation is a process, whose 
explanation is still absent. 

In contrast, the mixed phase model 0 — one of the 
most recent representatives of the second category — fo¬ 
cuses on the state of thermodynamical equilibrium. It 
constructs a non-ideal equation of state which includes 
part of the particle correlations and is able to show con¬ 
finement behavior. Although — in our opinion — this is 
a promising research direction, it fails to describe a very 
non-equilibrium hadronization, where the characteristic 
time of bound quark cluster formation is commensurate 
to the thermal equilibrium maintenance time. Also, di¬ 
quarks and baryons have not been included yet. 
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There are several chemical approaches to hadroniza- 
tion, showing an evolution time of a few to several fm/c. 
Assuming a massive, effectively non-relativistic quark 
matter, the hadronization can be described by utiliz¬ 
ing phenomenological confinement potentials or just the 
perturbative Coulomb-like potential. The elimination of 
color charge during this process is then due to medium 
effects, like string formation. The time evolution of the 
mixed quark-hadron system is described by a set of rate 
equations. Such an approach has been published a decade 
ago , and is recently under development in the frame¬ 
work of ALCOR model . 

Also in the framework of Nambu Jona-Lasinio (NJL) 
model hadronization cross sections have been obtained 
Since mesons, like pions and cr-mesons, are ele¬ 
mentary in this model, the confinement mechanism, an 
important detail of the microscopic process, remains un¬ 
explored. An alternative model, the chromodielectric 
model (CDM) originally proposed by Friedberg and Lee 
iQ, describes hadrons as clusters of quarks bound in 
a phenomenological scalar potential. While there have 
been a tantamount number of studies of stationary states 
of this model with quantum-mechanical quarks [19-22|, 
a real-time dynamical approach is quite recent [23-26|. 
The latter studies utilized semi-classical transport the¬ 
ory, treating quarks like classical charge clouds. 

In the present paper we report about a classical molec¬ 
ular dynamics approach on the quark level to hadroniza¬ 
tion. For the dynamical description, the classical limit 
of the chromodielectric model is used, where quarks and 
gluons are represented by extended but rigid charge dis¬ 
tributions. The background helds are a scalar held and 
two abelian helds, corresponding to the two neutral col¬ 
ors in SU{3). The six remaining, charged, gluons are 
treated as particles. We deal with these simplihcations 
in order to be able to carry out a molecular dynamical 
calculation of about 400 quarks in an acceptable time. 
The essential conhnement mechanism in this model is 
hidden in the scalar held self-interaction potential and 
the construction of the chromodielectric constant, whose 
particular dependence on the effective scalar held cr in¬ 
cludes the main quantum-mechanical effect from QCD. 

We think that such microscopical studies may re¬ 
veal whether or not the hadronization process is faster 
than the maintenance of thermal or chemical equilib¬ 
rium, whether or not the formed hadrons are close to the 
predictions of constituent quark models, and what kind 
of color-neutral clusters form in abundance. We think 
that the drawbacks of this approach, namely the clas¬ 
sical treatment of the color charge and its clusters, the 
abelian dominance assumption in the long-range fields, 
and the use of the cr-field, are outweighted by the advan¬ 
tages: an event by event, microscopic, dynamical simu¬ 
lation including flux-tubes that are non-elementary but 
formed by a dynamical process. This is useful for studies 
of possible hadronization scenarios of quark matter until 
real-time QCD calculations become feasible. 

This paper is organized as follows: in Sec. a brief 


description of the original chromodielectric model is fol¬ 
lowed by an explanation of the confinement mechanism 
in this model. Section 111 presents our practical approach 
to the model: we deal with the treatment of color charges 
in the classical limit, present our notion of hadrons within 
the model, and touch upon relevant computational meth¬ 
ods. In particular, the principles of a fast adaptive multi¬ 
grid finite element procedure for the solution of the Gauss 
law are discussed. Section IV reports on our choice of 
model parameters and its consequences on statical model 
properties. Section ^ is the main part of the present ar¬ 
ticle: it presents actual simulation runs of hadronization 
scenarios. We describe our choice of initial states and 
show snapshots of the program runs. We then compile 
our results on cluster-formation, and show mass, rapidity, 
and transverse momentum distributions for white clus¬ 
ters. Finally, Sec. VI summarizes our conclusions. 


II. THE CHROMODIELECTRIC MODEL 


Since free quarks have never been observed experimen¬ 
tally, it is generally accepted that any physical state of 
finite energy has to be a singlet under color-S'C/(3) trans¬ 
formations. This phenomenon is called color confine¬ 
ment. Since we know that classical (“tree-level”) chromo¬ 
dynamics does not enforce color confinement, this has to 
be a quantum (“loop”) effect, characteristic of a strongly 
interacting non-abelian gauge held theory. However, a 
proof of color conhnement based solely on the QCD La- 
grangian has not yet been given, mainly because it is too 
difficult to identify the relevant QCD degrees of freedom. 

The chromodielectric (Friedberg-Lee) model tries 
to mimic QCD’s color conhnement on the classical level. 
It introduces an additional scalar held a that enforces 
color conhnement in a conceivable way. 


A. Model Lagrangian 


The Lagrangian reads 

C = 

with 

( 

/ 

(2.3) 

= ^(9^cr)(a'"cr) - 17(cr). (2.4) 

The symbol q-^ denotes a quark spinor of havor /; both 
the color and Dirac indices are suppressed, and the Feyn- 


, + Cc-\- C-a ( 2 - 1 ) 

ip-mf)qf ( 2 . 2 ) 


2 










man dagger notation p = is used. The covari¬ 
ant derivative = id^ — contains the color vec¬ 

tor potential = A^k^-, where the A“ denote the eight 
Gell-Mann matrices. The color field tensor is defined as 

:= ^[D^,D,]. (2.5) 

^ 9v 

It obeys the classical chromodynamical field equation in 
a covariant medium with dielectric function k((t), 

= r, (2.6) 

with color current 

f = 

/ 

=:r“Y- (2.7) 

The quark field obeys a Dirac equation, 

(ip — mf)q^ = 0 ( 2 . 8 ) 

with a constituent quark mass m/. The quarks can in 
principle be treated quantum-mechanically here, while 
the color and a field are regarded as classical fields. 

The scalar field a obeys a Klein-Gordon equation with 
a fourth-order polynomial self-interaction U(a) and a 
complicated source term 

df,df^a + U'{a) = (2.9) 

Both functions U and k are drawn schematically in Fig. |^. 
Note that there is no direct coupling between the quarks 
and the tr—field in this model. 



FIG. 1. Schematical graph of the functions U{a) and K.{a). 
Note that U is an energy density while k is dimensionless. 
(Tvac denotes the value of the a field in the physical vacuum 
(ground state). 


B. Confinement Mechanism 


The Lagrangian (2.1) describes a classical field theory 
with an explicit color confinement mechanism. The in¬ 
gredients of this mechanism are the self-interaction term 
U{a), the dielectric function /c(cr), and the color Gauss 
law. In this section, we utilize a simple example to ex¬ 
plain how these three ingredients work together to en¬ 
force color confinement. 

First note that in the vacuum of this theory, a = cTyac 
holds everywhere and thus n = 0. 



Now consider the static and color neutral quark- 
antiquark configuration depicted in Fig. |^. The charge 
distribution is given by jo{x) = jo{x)^, i. e. it belongs 
to a U{1) subgroup of SU{3). The Gauss law reads for 
this case 

= jl (2.10) 

It forces a color flux = n{a) VAg — dtA^^ to 
stretch from the quark to the antiquark. This flux cannot 
exist in regions where k = 0, therefore the a field has to 
deviate from its vacuum value, forming a tube (or bag) 
with a < Uyac for the flux. Assuming U and k are de¬ 
fined so that cr ends up being zero in the bag, the bag 
formation process costs an energy density B, which in 
turn can be interpreted as the vacuum pressure exerted 
on the bag. 

Formulated differently, the vacuum pressure generated 
by U{a) and the coupling term kF'^ confines the color 
flux to a string — U and k together generate the color 
flux confinement. On the other hand, the quarks cannot 
escape from this string because the Gauss law demands 
a contiguous color flux to connect them. Unless pair cre¬ 
ation is allowed for, the resulting string can become ar¬ 
bitrarily long, causing the qq—potential to rise linearly. 

Let us contrast this with the corresponding result in 
classical chromodynamics. The Gauss law, reduced to 
the abelian subgroup of our problem, reads then 

V • (2.11) 

The solution is — like in electrodynamics — simply the 
linear superposition of the Goulomb fields of the two 
charges. Since the long-distance Coulomb potential is 
oc 1/r and not linear, we see that classical chromody¬ 
namics does indeed not include color confinement. 


3 















III. TREATMENT OF COLOR CHARGES, FIELD 
EQUATIONS, AND HADRONS 


A. A Classical Molecular Dynamics Approach for 
the Charges 


On our way to build a computer simulation of the chro¬ 
modielectric model, we have to make several approxima¬ 
tions. In doing so, we take care to preserve color confine¬ 
ment, since this is what we want to simulate. 

We treat only the two commuting gluon fields 

and Al^ as classical fields; they correspond to two 
decoupled, color neutral fields which behave just like 
Maxwell fields. So from now on, the color index a is 
generally restricted to the values 3 and 8. This approach 
amounts to taking a f7(l) x 17(1) abelian subgroup of the 
full SU{3) gauge group of QCD, and is therefore called 
the abelian approximation However, in contrast 

to the cited works, we do not simply drop the other, 
charged, gluons. Instead, they are treated like the quarks 
and antiquarks: as classical charge distributions with a 
Gaussian shape of fixed width. Spin and isospin are 
accounted for by appropriate deg ene racy factors. This 
changes the particle Lagrangian (2.2) into the classical 
expression 


^cjass ^ ( 3 . 1 ) 


with the color current 


f = gvY^PNix- (3.2) 


The index i runs over all quarks, antiquarks and charged 
gluons in the system. 


For the quarks, the color charges in (^) are equal 
to the diagonal entries in A^’®. Antiquarks carry the 
charge —q^'^, and the gluons carry each the sum of one 
quark and one antiquark color charge. Table | and Fig. ^ 
summarize the possible values. 

The masses mi are model parameters; our part icular 
choices for them are listed and explained in Sec. IVA . 
The function is a Gaussian number density distribu¬ 
tion of unit norm and RMS radius (x^) = 3rQ, 


PN{x-Xi) = 


27rr^ 


3/2 


,-{x-Xif/2rl 


(3.3) 


In the present version of our model, particle produc¬ 
tion or annihilation is not yet included. In principle, it is 
no essential problem to implement them, once an appro¬ 
priate probability density is at hand. This remains for 
future work. 



FIG. 3. Quark and gluon colors. The bold equilateral trian¬ 
gle corresponds to the quark colors, the other triangle to the 
antiquark colors, and the regular hexagon to the six charged 
gluons. The denotations r, g, b, f, g, b are chosen arbitrarily. 


TABLE I. Quark and gluon color charges in the abelian 
approximation. The color names are chosen arbitrarily but 
consistently throughout this article. The color shortcuts are 
the ones used in Fig. [^. 


particle 

color name 

q^ 

q^ 

q 

red r 

1 

1/V3 

q 

green g 

-1 

1/^3 

q 

blue b 

0 

-2/V3 

q 

antired f 

-1 

-l/\/3 

q 

antigreen g 

1 

-l/\/3 

q 

antiblue b 

0 

2/73 

9 

red-antigreen rg 

2 

0 

9 

red-antiblue rb 

1 

73 

9 

green-antiblue gb 

-1 

73 

9 

green-antired gf 

-2 

0 

9 

blue-antired bf 

-1 

73 

9 

blue-antigreen bg 

1 

73 
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B. Hadrons as Irreducible White Clusters 


If a quark-gluon plasma is formed in an ultrarelativis- 
tic heavy-ion collision, the droplet as a whole has to be a 
color-5't7(3) singlet. Since our numerical model respects 
only an abelian C/(l) x t/(l) subgroup of the SU{3) sym¬ 
metry group, an “SU (3) color singlet” is represented sim¬ 
ply by a “white” cluster, i. e. one that is color neutral 
with respect to both the A^— and the A®—charge (cf. Fig. 

An SU{S) color singlet is then built “automatically” 
by the presence of classical A® and A® fields, as solutions 
of ( |2.6D . Our simulation therefore starts off from a large 
but white quark-gluon plasma droplet. 

In nature, the hadronization process then divides the 
initial large cluster into a number of hadrons that are 
themselves separate color singlets. Correspondingly, our 
simulation is supposed to end up in a hnal state of small, 
white clusters that can be interpreted as hadrons. In 
this section, we investigate and classify the sorts of small 
clusters that we expect to find in the final state. 

Imagine a small, white cluster of color charges. Sup¬ 
pose it cannot be subdivided into smaller but still white 
subclusters. Then, flux tubes bind the cluster together, 
forcing it to rotate or oscillate in a “yoyo” mode 
Thus it will emit a waves, which cool the cluster down. 
Finally, all the color charges are resting with respect to 
each other, being centered at the same point in space. If 
we choose the same color charge distribution width for 
all particle species, the final cluster’s charge density will 
vanish identically; therefore, all helds will take on their 
repective vacuum values in this cluster. 

Because of the vanishing of all fields, two such clus¬ 
ters that happen to be close to each other can neither 
attract nor repel, and a small relative momentum suf¬ 
fices to separate them spatially. Consequently, we may 
conclude that the final state of our simulation must be a 
collection of irreducible white clusters, i. e. white clusters 
that cannot be divided into white subclusters. It turns 
out that only a hnite number of possible color combina¬ 
tions can form irreducible clusters, which enables us to 
classify them, according to Fig. ^ as mesonic, baryonic 
and pure glue ball states. The latter occur as mesonic 
resonance states in the experimental observation. 

We can turn this argument around, asking what clus¬ 
ters of quarks and gluons in our model correspond to 
hadrons. A hadron is supposed to be white but strongly 
bound. This leads us again to the notion of an irreducible 
white cluster. If the hadron is in its ground state, the 
colored constituents are all lying on top of each other, 
creating a local color charge density that is vanishing ev¬ 
erywhere. Consequently, a cold hadron in our model is 
completely field-free, and its mass equals simply the sum 
of the masses of the constituents, without any interaction 
contribution. E. g., look at the (qqq) baryon displayed in 
Fig. The three quarks carry the colors red, green, and 
blue, in the meaning of Fig. This baryon is certainly 
white, and it is also strongly bound: if one of the quarks. 


say the blue one (g® = 0, g® = —2/-\/3), tried to escape, 
a flux tube of a strong H®—field would build up, pulling 
the fragments together with a force of 1 GeV/fm. In 
the ground state, the three color charges sit right on top 
of each other, creating a locally vanishing color charge 
density. In this case, all helds take on their respective 
vacuum values, and the cold baryon possesses a mass 
3mq, where mq is the quark mass. 


irreducible mesonic clusters 
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irreducible baryonic clusters 



irreducible glue clusters 




FIG. 4. Classification of irreducible white clusters in the 
color diagram of figure ^ Points and circles indicate the oc¬ 
cupation number of each color state. Numbers in brackets 
indicate the color degeneracy of the diagram for each quark 
flavor combination. 


When studying hadronization processes, our program 
checks at each timestep whether there are irreducible 
white clusters in the system. A “cluster” in the computa¬ 
tional sense is a set of color charges located in a common 
a held bag (dehned here as a connected spatial domain 
with K > i). The b ag vanishes for cold clusters in our 
model (see Sec. Ill B| ) , so if a set of charges is outside any 
a held bag but very close in phase space, we treat it as 
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a cluster as well. If the cluster is irreducibly white, i. e. 
possesses one of the color structures displayed in Fig. 
it is considered to be a new hadron. The program then 
classifies the hadron according to its color structure, com¬ 
putes particle and field contributions to momentum and 
energy, removes the charges and field imprints from the 
simulation as well as possible, and generates a hadron, 
i. e. a white, non-interacting object that just proceeds 
on a straight path and may even leave the lattice. 

A short side remark is in order: In the investigation 
of the chromodielectric model, many authors (including 
the original inventors of the model) added a term 




f 


(3.4) 


in the Lagrangian, coupling the scalar quark density di¬ 
rectly to the cr field. This generates a high constituent 
quark mass outside a bag, which causes a high-density 
phase transition to occur in the constant field approxi¬ 
mation (see, e. g., p^ , p^ ). The coupling term also estab¬ 
lishes bound solitonic states with vanishing 

color charge density. Here we do not include the term 
(3.4) in the Lagrangian (2.1). If we did, the irreducible. 


final clusters of a simulation run would not be field-free 
but rather generate a local dip in the a field. This would 
generate a short-range attractive force between the white 
clusters: they would forever stick together, cool down by 
the emission of a waves, and end up as a single large and 
cold object, spoiling the hadronization process we want 
to model. 

In the model flavor presented here, there is neither an 
attractive nor repulsive classical force between cold (field- 
free) hadrons. It is interesting to note in this context 
that Koepf et al. have found a soft and short-ranged 
A^A^—repulsive potential in a quantum mechanical inves¬ 
tigation within certain approximations of the very same 
model. 


C. Numerical Solution of the Field Equations 


Let us recall the equations of motion obtained from the 
Lagrangian © with the classical particle term (0- 
The quarks and charged gluons feel a Lorentz force com¬ 
ing from the color field. 


mi- 


(Pxi 

W 


= Qz 


Ea 


dxi 

dt 


X 


(3.5) 


The a field equation, rewritten from (2.£), reads 


9^cr 

W 


= VV-C/'(a)- 


K'{a) 


- (.B“)2) . (3.6) 


The a field is stored on a 125^—sized lattice with a vol¬ 
ume of (25fm)^, and both the particle and cr field equa¬ 
tions of motion are solved by a simple explicit timestep 
method (staggered leapfrog, step-length 0.04 fm). 


The color fields obey the classical Maxwell equations, 
V • (3.7) 


V X = - 




V • H = 0 


V X + 


dK{a)E° 

Wt 


(3.8) 

(3.9) 

(3.10) 


with charge and current density given by (3.2). These 
equations cannot be treated by an explicit timestep 


method. The reason is that the Gauss law (3.7), being 
a constraint equation, contains no time derivative and 
therefore constitutes an implicit condition that has to be 
satisfied at each timestep separately. 

To simplify the problem, we neglect the fields in 

our simulation (and thereby neglect the corresponding 
color waves). Although breaking Lorentz invariance, it is 
not as bad an approximation as it may seem: it is exact 
for static charge configurations like in the computation of 
the string tension (Fig. |^; it is also exact for yoyo-like ex¬ 
citation modes of flux tubes in their rest frame . Fur- 
ther more, the color confinement mechanism explained in 
Sec. IIB is not broken by this approximation. Finally, 
the color field waves cannot mediate forces or transfer en¬ 
ergy between hadronic bags because the Poynting vector 
kE X B vanishes anyway in the physical vacuum. So the 
neglection of B amounts to a neglection of some inner 
degrees of freedom of clusters that are nonessential for 
the confinement transition. 


Neglecting i3“, it is clear from (3.8) that can now 
be derived from a single potential. 


= -V0“. 


(3.11) 


The Gauss law (3.7), rewritten in terms of (jf, reads 


V • (k((t)V()>“) = 


(3.12) 


This is still a constraint equation, modelling the con¬ 
straint of quark confinement. It constitutes a boundary 
value problem that has to be solved at each timestep of 
the simulation. This has to be done carefully since the 
color electric flux is an essential ingredient of the confine¬ 
ment mechanism. 


Note that (3.12) is singular in regions where n van¬ 
ishes; this is a major obstacle to numerical treatment. 
We overcome it in practice by ch oosing a no nzero but 
small value for n^ac (see Sec. IV A ). Still, (|3.12| ) is nearly 
singular and correspondingly difficult to solve. A well- 
tuned successive overrelaxation scheme (SOR) on a fixed 
grid with V® points, possessing at the very best a time 
complexity of 0{N'^) |3^, does not suffice here. Fourier 
transformation methods are excluded because of the non¬ 
constant coefficient function k(ct). So we choose a fast 
solution algorithm based on the hierarchical basis multi¬ 
grid method Q. We are using a tetrahedral mesh of 
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finite elements and take piecewise quadratic basis func¬ 
tions for the Ritz-Galerkin variational ansatz. During 
the multigrid iterations, the mesh is adaptively refined 
by bisection |^, which leads to block-diagonal matrices 
on each level of the multigrid scheme ||^]. Refinement is 
performed only in those regions where the solution func¬ 
tion varies with short wavelengths, e. g. close to the 
charges. This guarantees a certain precision of the solu¬ 
tion (jf' independently of the structure of the sources 
or the coefficient function k((t). 

There is a problem related to image charges. Impos¬ 
ing Dirichletian boundary conditions on the boundary 
surface of our cubic domain, any color charge that comes 
close to the boundary “sees” an image charge reflected at 
the boundary surface. In fact, even flux tubes between a 
charge and its image charge can be created, which yields 
a particularly spectacular picture if a single color charge 
is put in the model world. Image charges, however, are 
unphysical because the domain boundary is unphysical 
— real-world space is infinite. Fortunately, there is a 
simple trick to remove image charges: we make the cu¬ 
bic domain of the </)“ field much larger — (40fm)^ — 
than the (25fm)^ lattice of the cr—field. Since no color 
charge can leave the inner lattice (it would go into a re¬ 
gion where a = Oyac), it will never get close enough to 
the domain boundary to feel forces from its image charge. 
The adaptive mesh refinement method allows us to en¬ 
large the computational domain like this (a factor 5 in 
its volume!) without an even noticeable increase of com¬ 
puting time. 

Typical meshes resulting from the refinement process 
contain about simplices (tetrahedra) on about 20 

levels. The Ritz-Galerkin ansatz decays on the various 
levels into about lO'^ small linear systems of equations, 
each linear system involving about 40 nodal field values. 
This huge data structure has to be rebuilt from scratch in 
each timestep, since K{a) is changing. To save computing 
time, caching of matrices and lazy evaluation techniques 
prove to be extremely useful. 

The whole process of solving the Gauss law takes about 
30 steps, in which adaptive refinement alternates with re¬ 
laxation (solution) of the underlying equation. In each 
such relaxation run, the hierarchical multigrid method re¬ 
laxes the Ritz-Galerkin equations in some 10-30 V-cycle 
iterations; their convergence is further accelerated by an 
Aitken transformation method applied to all nodes^ _ 

The GPU time thus needed for one solution of ( 3.12| ) 
varies between 15s and 200s on a Pentium Pro 200MHz 
machine, depending mainly on the structure of the source 
function p® and At, which in turn depend mainly on the 
number of charges present. 

The quality of the solution is such that we never ob¬ 
serve free color charges in actual simulation runs, and 
even the interaction between larger clusters that carry 
nonzero total color is strong enough that they do not 
separate. 


IV. STATICAL MODEL PROPERTIES 

Before studying dynamical properties, like the time 
evolution of a quark-gluon plasma, in the present model, 
let us briefly discuss our choice of model parameters, 
which is based solely on certain desired time-independent 
observables, namely the hadronic mass spectrum, and 
color flux tube properties. 


A. Choice of Model Parameters 


We choose the same charge distribution radius 
= 0.7 fm for all the various particle flavors. The 
model possesses in this case a particularly simple “hadron 
spectrum”. Cold (unexcited) hadrons are white (color 
neutral) combinations of several particles centered at the 
same location. Their color charge distribution vanishes 
everywhere, therefore the fields all take on their vacuum 
values. Thus, the mass of such a hadron is simply the 
sum of the masses of its constituents. In order to obtain 
roughly the proper hadronic mass scales, we have chosen 
a nonstrange quark mass of 400 MeV, yielding a mass of 
800 MeV for the p—meson and 1200 MeV for the nucleon. 
The other masses are all given in table |^. This table 
exhibits a reasonable agreement between model and ex¬ 
perimental numbers; fine-tuning could certainly improve 
this even further. 

TABLE II. Quark and gluon effective masses, and corre¬ 
sponding hadron masses [Particle Data Table], All masses 
are measured in MeV. The symbol q denotes an up or down 
quark, s a strange quark, and g a gluon. 


particle 

model mass 

exp. mass 

error 

q 

400 

- 

- 

s 

550 

- 

- 

g 

700 

- 

- 

i(A-f A) [qqq) 

1200 

1086 

-flO.5% 

Y [qqs) 

1350 

1385 

-2.5% 

H (qss) 

1500 

1530 

-2.0% 

D (sss) 

1650 

1672 

-1.3% 

p {qq) 

800 

770 

+3.9% 

LO (qq) 

800 

782 

+2.3% 

K* {qs) 

950 

892 

+6.5% 

>!> (ss) 

1100 

1020 

+7.8% 

99 

1400 

1400 ? 1 

- 
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The parameters of the quartic potential U (a) and the 
dielectric function K{a) are harder to fix since there is no 
well-known physical particle associated with a. However, 
if we interpret ct— waves as pure glueballs with a mass 
of about IdOOMey, we obtain U"{ayac) = (IdOOMeH)^. 
Furthermore, the bag constant B = {I'bOMeV)^, and 
CTjjac = 0.31fm“^ are chosen in order to get a low bag 
surface tension — which is easier to handle numerically 
— while still having a region with U"{a) < 0, which is a 
necessary condition for a first order color deconfinement 
phase transition to occur. Our particular choices are mo¬ 
tivated by studies of solitonic states in the model . 

Finally, we use 

U{a) = B + aa^ + ba^ + ca'^ (4.1) 

with a — 6.184fm~^, b = —80.72fm~^, c = 163.1, 
B = (150MeV)"‘. For k, it turns out that the simula¬ 
tion results are not very sensitive to its precise shape. 
As described in the last section, its only purpose is to 
couple the vacuum pressure generated by U{a) to the 
color flux tube. For doing this, the k function has to 
take on a value much smaller than one near a « <Jyac- 
Secondly, k has to rise quickly to one where U{a) ~ B, 
supporting our interpretation of B as bag constant. We 
use 

K(cr) = -- (4.2) 




-f 1 


with a = 5 and f3 = 0.4. 

Finally, the strong coupling constant as = is 

chosen such that the qq string tension equals about 1 
GeV/fm, a value that can be deduced from quarkonium 
spectroscopy |^. Anticipating the results of the next 
section, we show in Fig. g that a value of as = 2 results 
in a satisfactory numerical string tension. 



length [fm] 

FIG. 5. The qq-flux tube energy (a and color field contri¬ 
butions) plotted against the separation distance of the quarks. 
The slope is about 1 GeV/fm. 


B. Statical Properties of Flux Tubes 


Of course, the string tension, plotted in Fig. || and al¬ 
ready discussed in the last section, is the flux tube prop¬ 
erty most relevant to the hadronization process. Also, all 
our model parameters are already fixed by the consider¬ 
ations done in the last section. Nevertheless, we would 
like to show pictures of color flux tubes and their fields. 
We do this in Fig. ||: Picture (a) shows the color flux 
£>8 = —K(a)V(/>^ between a blue quark and an antiblue 
antiquark (in the nomenclature of Fig. The vector 
field is displayed in a plane containing both quark cen¬ 
ters — the computation has of course been performed in 
three dimensions. The spheres at both ends of the flux 
tube represent the quarks; their radius indicates the RMS 
radius of the color charge distributions (0.7 fm). It can be 
clearly seen that the color flux is indeed nicely parallel, 
as would be expected from a string conhguration. The 
string radius turns out to be rather large — about 1 fm. 
In terms of our irreducible white cluster classification of 
Fig. 1^, this string is an excited meson of the gg—type. 

Picture (b) of Fig. ^ shows the a field for the same 
charge configuration. The contour lines are drawn at ct 
values of 0, 0.1, 0.2, 0.3 (from the inside outwards), where 
o'vac ~ 0.31. Remember that the only source term of the 
CT field is given by k'(ct)F'^/ 4; thus, the vector field drawn 
in picture (a) effects a field imprint that does not fully 
extend over th e qu arks. Now we added a direct quark-tr 
coupling term (^^), taking qs = 8 and approximating qq 
by the quark number density given in (^^). The result 
is picture (c), differing from (b) only slightly, since the 
quark-CT coupling affects only the endpoints of the flrrx 
tube. 

Lastly, picture (d) shows the ct field for a gluonic string. 
The total color flux is larger here: since the gluons carry 
each the sum of a quark and an antiquark color, the total 
gluon charge \/{q^)^~+Jq^ is larger by a factor of -v/S in 
comparison with a quark charge (see Sec. Ill A ; for this 
picture, we took the gluon colors rq and qf). The total 
color fliDc is therefore enhanced by the same factor (en¬ 
forced by the Gauss law). This results in a larger pressure 
of the flux lines, which dig a significantly wider “hole” in 
the nonperturbative vacuum. In terms of our irreducible 
cluster classification scheme of Fig. |[ this gluonic string 
is an excited glueball of type qq. 


2 



-2 




















FIG. 6. (a): Vector field of a quark string of length 4 fm, 
(b) its a field, (c) modification caused by an additional scalar 
coupling term (3.4) with gs = 8, (d) a field of a gluonic string. 
The contour lines are defined by a = 0, 0.1, 0.2, and 0.3 
{(Tvac ~ 0.31). Axis ticks are in units of 1 fm. 


We can also use our program to investigate a possible 
string-string interaction. The cr field of parallel and an¬ 
tiparallel configurations of two strings is plotted in Figs. 
and For the parallel case (Fig. |^), we see no notable 
deformation of the strings as they ^proach each other, 
while in the antiparallel case (Fig. the strings “flip”, 
since the field configuration of lowest energy (which is 
the static solution of the field equations) is the one that 
connects only the closest quarks by flux tubes. 

Figs. H and 10 show the corresponding string-string 
potentials. For the parallel case (Fig. the potential 
differs little from a constant, possessing a slightly attrac¬ 
tive part (bag fusion) and a repulsive core (Coulomb re- 
pulsion). The potential of the antiparallel configuration 
(Fig. 1^ shows the expected stringlike (linear) behaviour 
as soon as the strings are flipped (below 4 fm). 





T 




FIG. 7. Fusion of parallel flux tubes. The contour lines are 
at the same a field values as in Fig. |^. 
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FIG. 9. Flux tubes flip in an antiparallel configuration. 
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FIG. 10. String-string potential for the antiparallel case. 


V. HADRONIZATION OF A QUARK-GLUON 
PLASMA 

The described simulation tool enables us to follow the 
microscopic dynamics of hadronization. However, in or¬ 
der to gain a qualitative insight, we restrict the first 
analysis to a few observables. In the next section, we 
present our choice of hadronization scenarios. Then, we 
introduce some important observables and explain our 
means of obtaining them from simulation data. Finally, 
we discuss our results on mass distributions, the evolu¬ 
tion of the quark-gluon-hadron matter composition, and 
analyze spectral properties. 


A. The Initial State of the Simulation 

When carrying out a molecular dynamic simulation 
on the quark-gluon level, one has to make compromises. 
Even using modern computers and state-of-the-art nu¬ 
merical methods, we cannot simulate systems larger than 
a few hundred of particles in a volume of (25fm)^. The 
entire simulation run takes finally about a 200 fm, cor¬ 
responding to weeks of CPU time on a Pentium Pro 200 
MHz system. 

Due to the above mentioned restrictions we decided 
to simulate only a central part of the heavy-ion collision 
region, including about 400 quarks and gluons in an ini¬ 
tial volume of about 100 fm^. In real Pb Pb collisons 
at CERN SPS, about 3000 — 4000 particles are present 
in the quark matter phase. We compare two different 
scenarios in order to gain a qualitative insight about the 
hadronization: (i) a fully stopped quark matter at ini¬ 
tial temperature Tq = 160 MeV, and (ii) a hot quark- 
gluon plasma initialized on top of a Bjorken-flow stretch¬ 
ing from coordinate-rapidity r/o = —0.5 to ijq = 4-0.5. 

The reaction volume is calculated assuming a radius of 
the initial sphere of i?o = 4 fm or of the initial cylinder 
of i?o = 4 fm, respectively. The initial energy density 
of Eo = 2.5 GeV/fm^ is obtained from the total bom¬ 
barding energy characteristic for SPS. The initial num¬ 
bers of quarks and gluons are chosen according to non- 
relativistic Maxwell-Boltzmann statistics, using the re¬ 
spective masses and the initial temperature. The ratios 
of various flavor quarks and gluons are determined from 
equilibrium principles (law of mass action), being pro¬ 
portional to 

/ T'’ \ 3/2 

iV,(0) oc , (5.1) 

with di the spin, color and flavor degeneracy. The re¬ 
sulting parton numbers are 133 light quarks, 40 strange 
quarks, the same numbers of respective antiquarks, and 
44 gluons for the Bjorken event, and 196 light quarks, 52 
strange quarks, and 57 gluons for the full stopping event. 
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The random Gaussian momenta (on top of the initialized 
flow pattern) satisfy 


^ = To (5.2) 

for each particle species and momentum components fi = 
x,y or z. We do not consider a nonzero quark chemical 
potential in the present work; furthermore there is no net 
color and no net strangeness in the system. 

A short side remark: A possible third scenario, a fully 
“transparent” picture where the particles are initially 
separated in phase space, leads to a rapid separation 
into unconnected parts (the spectators), leaving only a 
small fraction of mini-clusters in the mid region. Unfor¬ 
tunately, we cannot simulate too large relative momenta 
at present, due to the computational restriction to inter¬ 
mediate size simulation volumes. Some particles bounce 
back from the lattice boundary and interfere artificially 
with others still in the expanding fireball. This is unphys¬ 
ical, and we have to abandon our studies of this scenario 
at this point. 



B. Identification of Observables 


The basis of our identification of observables is the 
white cluster classification and hadronization mechanism 
described in Sec. IIIB. The total momentum and total 


energy of white clusters can be calculated from the in¬ 
dividual four-momenta of its constituents and from the 
field energy carried by the classical vector fields. The 
mass of a cluster is defined by 


-^clus “ ( T/field 


- 


Pfield 


E 


Pi 


(5.3) 


Here, 

T/field — 
Pfield — 


d^x ( -b ^(Vcr)^ -b U{a) + 


= — d'^x&Va. 


(5.4) 


The mass distribution, dN/dM, tells the number of 
clusters with a mass between M and M -b dM. The reso¬ 
lution dM of the mass binning is chosen to be dM = 100 
MeV presently. If dM is too large, one cannot distin¬ 
guish different particles; if it is too small, the statistics 
becomes too low. The above value is a compromise found 
by trial and error. 

Further observables are the rapidity distribution 
dN/dy with 



E + Pz 

E-Pz' 


(5.5) 


and the transverse momentum distribution dN/dpr with 


PT = \JpI+pI- (5-6) 

Collective flow and local temperature both affect these 
spectra: a strong longitudinal flow causes a plateau in 
the dN/dy distribution, while the slope of the transverse 
momentum spectra — if nearly exponential — includes 
both thermal and radially collective motion. 

Finally, some global variables are of interest in order to 
obtain an approximate equation of state from the molec¬ 
ular dynamic simulation. A temperature-like quantity 
can be obtained by fitting the non-relativistic Maxwell- 
Boltzmann distribution to the momenta by inspecting 

= {pl^)/m^, (5.7) 

for a cluster of type i with mass rui, and momentum di¬ 
rection p = x,y or z. The approximate equality of these 
temperature-like parameters for different spatial direc¬ 
tions p reflects isotropy of the pressure, while its equality 
between different sorts of clusters i chemical and thermal 
equilbrium. 

It remains for future work to analyze the cluster size 
distribution with its higher moments after collecting im¬ 
proved statistics. This can shed some light onto the order 
of the hadronization transition from the viewpoint of per¬ 
colation theory, as recently pointed out by H. Satz Q. 
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FIG. 11. Simulation snapshots after initialization of a central-rapidity slice of a QGP tube with Lorentz-invariant flow 
(“Bjorken scenario”). The frame diagonal is of length 45 fm. The labels at the particles and hadrons indicate the various 
flavors and constituents; the letter q represents an up or down quark; the other letters have obvious meanings, capital letters 
denoting antiquarks. Colored objects (quarks, gluons) are drawn in dark grey, hadrons in light grey. 
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FIG. 12. Simulation snapshots after initialization of a thermal QGP fireball (“full stopping scenario”). 
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C. Example Simulation Snapshots 


Figs. 11 and show six snapshots of the time evolu¬ 
tion for the Bjorken and “full stopping” scenario^ The 
simulation times of the snapshots are written inside the 
frames in units of fm; the frame diagonal is of length 45 
fm. Quarks, gluons, and hadrons are drawn as spheres; 
for better visibility of the scene, the radius of a colored 
object is scaled down to one third of the actual RMS 
radius of the charge distribution (0.7 fm). Hadron radii 
reflect the bag volume formed by the a field footprint 
of the hadronized cluster; they are scaled down to one 
half of their actual size. For white clusters without any 
a field bag (see section HI B), we draw a sphere of radius 
0.35 fm. The a field as well as the two color fields are 
not displayed. Note that in the last frame, many hadrons 
have already left the scene. 

As can be seen from the pictures, the initial partonic 
fireball is expanding and hadronizing. At each timestep, 
the equations of motion of the charges and t he a field as 
well as the Gauss law is solved (see Sec. |HIC ). Whenever 
an irreducible white c luster is detected (by the mecha¬ 
nism described in Sec. HIB), it is classified according to 
Fig. H as a meson, baryon or glueball, its total energy, mo¬ 
mentum and mass are computed, the partons and their 
field imprint are removed from the simulation and re¬ 
placed by a noninteracting “hadron”, displayed in the 
snapshots as a grey sphere. We would like to emphasize 
that the production of hadronic states is not artificially 
put in by hand: the model itself produces irreducible 
white clusters in the final state, and it is only up to our 
detection routine to classify them and take them out of 
the simulation. The detection routine does not drive or 
interfere with the process of hadron formation, it merely 
serves to interpret the results. Although the last snap¬ 
shot still shows quarks and gluons, this is not the final 
state of the simulation run. In the final state of the simu¬ 
lation run, hadronization is indeed complete: no partonic 
fireball is left over. 


D. Energy Conservation, Hadronization Time, and 
Temperature 


Fig. 13 shows the energy over time of the system and 
its parts in the Bjorken scenario, for the same event the 
snapshots were taken from. Fig. 0 shows for the same 
time interval the temperature evolution of the partons in 


^By the time this article will be printed, MPEG movies 
of the actual simul ation ru ns will be retrievable by anony¬ 
mous ftp from ftp://theorie.physik.uni-giessen.de 


/pub/cdmmovie/bjorken.mpg resp. /stopping.mpg 


the fireball and the produced hadrons. Their tempera¬ 
ture is found from 



The partons heat up initially because they are accel¬ 
erated by relatively long (> 1 fm) flux tubes. These flux 
tubes are present simply because of the initial spatial 
random distribution of the charges. Very soon, however, 
the hot particles cool down, transferring energy to the a 
field. On the same time scale, the energy contribution of 
the color electric field, also shown in Fig. |^, decreases to 
almost zero. This is an indication of the strong electric 
shielding of the plasma. The picture of the plasma as 
a collection of flux tubes breaks down here; flux tubes 
emerge only on the surface of the fireball, pulling back 
colored objects that try to escape. It is interesting to note 
that, although the energy contribution of the color field 
itself is never significant, the color field dictates the dy¬ 
namics of the hadronization process, allowing only white 
objects to separate. 

The hadronic energy contribution is a monotonous step 
function: it increases each time we take another quark- 
gluon cluster, delete its field imprint and add it to the 
set of hadrons in the system. The hadrons do no longer 
interact, hence the piecewise constant energy. About half 
of the final energy is deposited in a waves. This energy 
has its origin in the initialized thermal energy of the color 
charges. These form flux tubes, and via the coupling term 
«:(cr)F^/4 the flux tube cools down, radiating a waves. 
The (T—field describes the long-range collective effects of 
QCD. Its large energy content shows that the nonper- 
turbative QCD vacuum is highly excited; in the “real 
world”, it would decay mostly into pions which are so far 
not contained in our simulation. 

Note that most of the final hadronic energy is produced 
in the time interval [lOfm, 50fm]. These 40 fm could be 
interpreted as the hadronization time in our model. 

There is no final-state interaction among the hadrons 
or between hadrons and partons in the simulation, so 
the hadronic subsystem is not in thermal contact with 
the fireball. This is realistic: the hadron temperature 
reflects the parton temperature at the time of hadron 
formation (freeze-out point). Since most of the hadrons 
are produced in the time interval [lOfm, 50fm], the final 
hadronic temperature reflects the fireball temperature at 
that time. The fireball cools further down, which leads to 
the unrealistic picture that after about 50 fm, the fireball 
is significantly cooler than the hadrons. 

In an experiment, a fireball as cold as 50MeV would 
probably hadronize immediately. In our simulation, how¬ 
ever, hadrons are only well-defined once they have clearly 
left the fireball and can be detected as a separate, irre¬ 
ducible white cluster. At this late stage of the simulation 
run, hadronization is a random evaporation process: in¬ 
side the fireball, color fields are strongly shielded, which 
prevents stable white clusters to form. So an irreducible 
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white cluster (hadron) can only form at the surface of 
the fireball, as soon as its constituents, floating around 
in the fireball nearly freely, happen to reach the surface 
with coalescent momenta. Thus it takes quite some time 
for a small, cold parton cluster to hadronize completely. 

We attribute this partly to the missing repulive inter¬ 
action between the white clusters in our simulation 
In addition, it is well known from nuclear fragmentation 
studies that the mean field tends to prevent clusteriza¬ 
tion 1^. Thus, fluctuations of the color fields may be 
important for the final hadronization process. 



0 - --- ---^ 

0 100 200 


time [fin] 

FIG. 13. System energy during time evolution in the 
Bjorken event. 



0 ^^^^ 

0 100 200 


time [fin[ 

FIG. 14. Particle and hadron temperature over time in the 
Bjorken event. 


E. Mass distribution 

The final mass distribution of the hadrons is shown 
in Fig. 1^ for both the Bjorken and the “full stopping” 
scenario. The spectra start at a mass of 800 MeV because 
this is the mass of the lightest hadron [qq) in our model, 
and no quarks or gluons are left over in the final state. 

Although the spectra do look different on first sight, 
one should keep in mind that since these are single events, 
the differences in the spectra may simply be due to sta¬ 
tistical fluctuations and might not be significant of the 
respective scenario. Note also that the shape of the mass 
spectrum in the Bjorken case would stay the same if the 
chosen rapidity slice Ay = [—0.5, 0.5] were extended to 
a realistic interval, simply because masses are Lorentz 
invariant. 



0 12 3 4 


mass [GeV[ 



FIG. 15. Mass spectra; the upper graph is the Bjorken case, 
the lower the “full stopping” scenario. 
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Assuming a Hagedorn density of hadronic states , 


F. Rapidity and Momentum Spectra 


p{m) (X (5.9) 

we expect in the nonrelativistic limit a measured mass 
distribution 


^ oc p{m) f 
dm J 

oc (5.10) 

The nonrelativistic limit is justified by our large con¬ 
stituent masses. Here, Th is the Hagedorn temperature: 
the limit temperature of hadronic strings. Beyond this 
temperature, the flux tube picture must be replaced by 
a true QGP. We assume Th = 160 MeV for now (which 
is what Hagedorn estimated originally) and take a = —3 
f rom [^ . Fig. [Tsl shows that a satisfactory fit of the form 
( ^.10 ) can be obtained. Our fit parameter b is related to 
T via 


T Th' 


(5.11) 


The corresponding hadron temperature is T = 130 MeV 
for the Bjorken event (T = 146 MeV for the full stop¬ 
ping event), which is somewhat higher than the kinetic 
temperature found in the last section. However, these 
numbers depend on the choice of a: taking a = —3/2, we 
expect a simple exponential curve for the mass spectrum, 
which results in different values for T, namely T = 110 
MeV for the Bjorken event and 125 MeV for the full stop¬ 
ping event. In this single-event analysis with our very 
small statistics, we cannot really decide on the correct 
value of a. 

Another possible expectation for the mass spectrum 
could be that of a percolative process near the critical 
point iQ. There, a scaling distribution 


dN 

—— oc m 
dm 


(5.12) 




-1 0 1 

y 


is ex pec ted. Trying to fit this relation to our data (see 
Fig. 1^), we find a satisfactory agreement for r = 3.06 
(Bjorken case) and r = 2.22 (full stopping event). We 
plan to investigate this further in the future. 


FIG. 16. Rapidity spectra; the upper graph is the Bjorken 
case, the lower the “full stopping” scenario. 

Fig. ^ shows the hadron rapidity spectra of both 
scenarios. They do not look significantly different, the 
Bjorken case being somewhat broader due to the longi¬ 
tudinal flow. However, no plateau structure is visible in 
the Bjorken case, since we have cut out a central-rapidity 
slice of the Bjorken tube. One could easily extrapolate 
this result by folding the rapidity spectrum with that of 
a Bjorken QGP tube of realistic extension. This folding 
process would then result in a plateau structure. 

Finally, Fig. 0 shows the transverse momentum spec¬ 
tra of nonstrange white hadronic clusters in the Bjorken 
scenario, observed in mass bins of width 0.6 GeV above 
0.6 GeV. Although the points in the Figure representing 
mid-pr values do not lie exactly on exponential curves. 
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it is interesting to consider the slope values in order to 
extract possible information on the flow and temperature 
of the fireball in the late hadronization phase. The re¬ 
sulting inverse slopes, denoted as “temperatures” T{m), 
are displayed in Fig. Ill- By fitting to the linear relation 

T (m) = To + m{v^) ^ (5.13) 

we extract from the plot a temperature Tq = 77 MeV 
and a radial flow \/(?^ = 0.32. A direct compu¬ 
tation of the radial flow from our final state yields 
y/(v^) = 0.30, which is remarkably close. For the full 
stopping event, the same T(m) analysis yields To = 126 
MeV and y/{v"^) = 0.18, while the true radial flow is 
y/(v^) = 0.27, a not so good agreement. Compared 
to data on SPS heavy ion experiments at CERN 
{{vt) = 0.4, Umax = 0.6 and T = 140 MeV) our results 
display a little less flow and a colder hadronic fireball. 
It should be noted, however, that this analysis is based 
on a single event whose plasma volume is only a fraction 
of the experimental one. Also, we extracted the spectra 
of rather heavy hadronic clusters; their decay products 
(with lower mass) may have higher p^-slopes leading al¬ 
together to a higher estimate of Tq- From this viewpoint, 
the qualitative agreement is satisfactory. 



FIG. 17. Final transversal momentum speetrum for the 
Bjorken case, for three mass bins separately. When comput¬ 
ing the fit lines, we took the statistical error bars correctly into 
account 



FIG. 18. The temperatures extracted from mass differential 
Pt spectra, plotted against hadron rest mass. 


VI. CONCLUSION 

We have investigated the hadronization process of a 
quark-gluon droplet in the framework of the chromodi¬ 
electric model, using a classical molecular dynamic ap¬ 
proach. An efficient computer code was written and put 
to work for the solution of the boundary value problem 
connected to the Gauss law in the abelian dominance 
approximation. This code made it possible to study a 
volume of about (25fm)^, filled with about 400 particles. 

Our quark masses were fitted to a constituent model 
of the vector meson octet and baryon decuplet states. 
We determined from the microscopic data of hadroniza¬ 
tion events the mass-, rapidity- and transverse momen¬ 
tum distributions. 

As example scenarios we have taken two different ini¬ 
tial states, one with full stopping and one with a one¬ 
dimensional scaling flow. The two scenarios behave sim¬ 
ilarly, having a commensurate energy density initially. 
We found that although the identification of white clus¬ 
ters technically requires a long simulation time, the mass 
distribution develops a hadron part already at the begin¬ 
ning of the process — about 10 fm after initialization — 
and the main hadronization process is done after another 
25 fm have passed. In a recent RQMD paper |Q| similar 
hadronization times are reported. 

Hadronization is certainly a much faster process than 
this simulation suggests. This is because we can only 
uniquely identify an irreducible white cluster after it has 
already left the fireball. So in this simulation, hadroniza¬ 
tion looks like a slow evaporation effect, while in princi¬ 
ple pre-hadrons could already form in the plasma but re- 
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main undetected by our program. An effective, quantum- 
mechanical, repulsive force between white clusters, as 
computed by Koepf et al. |^|, could change this picture, 
if a consistent way of including this in our simulation can 
be found. 

Summarizing, we do not observe a vast departure from 
thermal equilibrium during this microscopic study of 
hadronization. Chemical equilibration on the other hand 
is incomplete as shown by the deviation between the 
Maxwell-Boltzmann slopes for different cluster species 
and in general between colored and white particles. 
This also supports the conclusion that during the late 


hadi'uiiizHLiuii of quark iiiattei' Lliere is no phase equilib¬ 


rium between the quaik-gluuii and hadiuii cumpuiieiiLs, 
the Gibbs criteria do not apply. The suggested picture 
is rather a mixture of hadronic and quark-gluon plasma 
counterparts (with additive instead of equal pressure). 

An extended study of cluster distribution in the near 
future is planned in order to investigate whether the 
hadronization is a percolation- or phase-transition-like 
process. 
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